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A QCD-based treatment of projectile size fluctuations is used to compute in- 
elastic diffractive cross sections (Jdiff for coherent hadron- nuclear processes. 
We find that fluctuations near the average size give the major contribution to 
the cross section with < few% contribution from small size configurations. 
The computed values of cr^iff are consistent with the limited available data. 
The importance of coherent diffraction studies for a wide range of projec- 
tiles for high energy Fermilab fixed target experiments is emphasized. The 
implications of these significant color fluctuations for relativistic heavy ion 
collisions are discussed. 
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We demonstrate that color fluctuations in the projectile wave function play an impor- 
tant role in high-energy (Eh) hadron (h)-nuclear reactions. This is done by studying the 
cross section Udiff for inelastic coherent nuclear diffraction which would vanish without such 
fluctuations. 

Not much data for inelastic coherent nuclear diffraction exists. But making the necessary 
measurements seems much more feasible now than 10-15 years ago, due to the development 
of micro-vertex detectors |I[ in which the target is the detector itself. Using these detectors 
can insure that the nucleus is left in its ground state in the diffractive process. Furthermore, 
measuring the coherent nucleon-nuclear diffractive cross section will ultimately allow the 
computation of fluctuations of various observables available in heavy ion collisions [0,0. 
This means that new experiments, which can be done during the 1995 Fermilab fixed target 
run, have gained a new urgency. 

We begin our analysis with a general discussion. It has long been known that the time 
scale, given by the uncertainty principle, for quantum fluctuations from a hadronic state h 
into a state X is % Eh 3 . Such fluctuations are inhibited at large enough energies, so that 
one may treat the hadron as frozen in its initial configuration [||. The natural approach to 
describe such collisions is the scattering-eigenstates formalism, introduced by Feinberg and 
Pomeranchuk || and Good and Walker ||. This accounts for the high energy coherence 
effects: the projectile can be treated as a coherent superposition of scattering eigenstates, 
each with an eigenvalue a. The probability that a given configuration interacts with a 
nucleon with a total cross section a is P(o). It is possible to reinterpret P(o) by relating 
the size of a given configuration with its cross section (forward scattering amplitude) in a 
monotonic fashion; see e.g. |T|-|TU| and the references therein. 

The scattering-eigenstates formalism accounts for the effect that different components of 
a hadron have different interactions with the target. For example: (1) In color transparency 
physics a configuration with a small size has small interactions. (2) Nucleonic configurations 
with and without a pionic cloud have different interactions with a target. (3) In string models 
of hadrons the transverse size is approximately proportional to the sine of the angle between 
the momentum and string directions. 

In quantum mechanics a system fluctuates amongst its components. Such fluctuations 
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lead to different kinds of interactions. Since color dynamics determines all the interactions, 
we use the term color fluctuations. 

What is known about the -P(cr)? The convenient procedure is to consider moments of this 
quantity: (a n ) = J daa n P(a). The zero'th moment is unity, by conservation of probability, 
and the first is the total cross section. The analysis of diffractive dissociation data from 



the nucleon 111 as well as inelastic corrections to the total hadron-deuteron cross section 



T2|j determined the second moment of P(cr) (see the summary in ||T3[) while diffractive 
dissociation data from the deuteron targets determined the third moment of -P(cr) for protons 
13| . The functional form of -P(cx) was then determined by taking the behavior for small 



values of a from QCD fll4| , |l5l and also including the rapid decrease of -P(cr) for large values 



of a. Thus it is possible to obtain a realistic form of P(cr) for a wide range of a [|T^,[nj. 

The ideas behind the formulae for adiff(A) were suggested a long time ago, see e.g. Refs. 
[ l~T|,p~6| p~8j| and Ref. [|19| and references therein. Ref. [20] reviews the attempts to describe 



data in terms of a few scattering eigenstates. However, a realistic model, based on QCD, for 
the cross section fluctuations was missing. For example, preQCD models contained terms 
corresponding to a 5(a) piece of -P(cr) [0,0 which QCD does not allow. Here we provide 
an expression for adiff(A) in terms of the independent information given by P{<j). We start 
with the standard formula for the diffractive cross section in terms of the transition matrix 
f: 

da diff (A) oc J2 8 \ P f -Pi)---\(A;a,M%\f\ A, h) | 2 , (1) 

where the ellipsis represents the standard phase space and flux factors. The frozen approx- 
imation allows us to use completeness to sum over the diffractive excitations a^M\. Then 
the only hadronic information resides in the square of the hadronic wave function. The key 
step, introduced by Miettinen and Pumplin [JTTJ] and revived in Refs. p|, p!3|JT5[| is to re-express 



the integral over that squared wave function in terms of an integral over a that involves the 
independently determined probability P(cr). For coherent nuclear processes the scattering 
wave function can be obtained using the optical potential, now also a function of the inte- 
gration variable a. We consider values of A greater than about 10, so that the t-dependence 
of the nuclear form factor is much more important than that of the hN diffractive amplitude. 
Then the coherent nuclear diffractive cross section adiff(A) can be expressed as: 



a diff (A) = I d 2 B f daP(a) £ [< h \ F(a, B) \ n > 2 ] - f daP(a) < h \ F(cr, B)\h> 

(2) 

where F(a, B) = 1 — e~* oT ( B } and T(B) = /f^ pa(B, Z)dz. Here the direction of the beam 
is Z and the distance between the projectile and the nuclear center is R = B + ZZ. Eq.(2) 



is similar to the one used in [19| which did not introduce the notion of P(<j). The advantage 
of eq. (2) as compared to the related equation of Ref. [^T[] (for a review see |2"2]] ) within 
the two gluon exchange model is that we do not need to assume the validity of pQCD at 
average interquark distances in hadrons where as is large (with an extra prescription for 
dealing with gauge noninvariant effects due to introduction of nonzero gluon mass), nor use 
constituent quark model wave functions instead of parton wave functions. 

It is instructive to consider the extreme black disk (bd) limit of eq. (2). Then the function 
F(a, B) is unity for positions inside the nucleus and zero otherwise, so that a d iff(A) vanishes! 
In particular, the black disk model gives = 2ttR\, a b e f = irR\ and a^jf(A) = 0. But we 
shall see that including the effects of color fluctuations leads to observable diffractive cross 
sections rapidly increasing with A, which are consistent with existing measurements of semi- 
inclusive diffraction [p3|24 |. Another way to show that color fluctuations cause adiff(A) is 
to observe that taking P(a) to be a delta function, e.g. P(cr) = S(a — a) gives ad%ff(A) = 0. 

We first display results for the pion projectile and use three parametrizations of P(cr) 

(<?-<?o) r - 



of Blattel et al |l5j of the form P(cr) = N(a,n)e ( ntT o) n . All of these distributions have 
approximately the same value of 



u n = 



(ay 



(3) 



with LU a = 0.4 — 0.5. The resulting cr di ff(A) are shown in Fig. la. Note that for each P(cr), 
adiff(A) varies as A im for A »s 16 and as A 0S5 for large A ps 200. 

Next we examine the adiff(A) that can be observed in proton scattering. The current 



data indicate |13| that (a ) w 1.25(a) (for proton energies of about 400 GeV) so we may 



expect interesting effects. The results are shown in Fig. lb for three versions of P(o) of 



ref. [13]. The parameterization for the proton is P(cr) = N(a, . The curves 

are obtained with n=2, n=6 and n=10. We see that the shape of P(cr) plays an important 
role in obtaining the magnitude of <Tdiff(A). The A-dependence is of the approximate form 



(<^o) r 



A 80 for A ~ 16 and A 4 for A ~ 200, which is smaller than for the pion case because here 
P(a = 0) = and the average value of a is larger. 

At present there are no data available on the A-dependence of a diff {A). However the 
A-dependence of the reaction 7r + + A — ► n + + n + + 7r _ + A was studied in p3[ for p n + = 200 
GeV. This data integrated over the mass interval 0.8 < < l.bGeV (and corrected for 
the small contribution of Coulomb excitations) are shown in Fig. la. In Fig. lb we present 



the data of |24|]on n + A — ► pn~ + A for the the mass interval 1.35 < M < lA5GeV 
with the Coulomb contribution subtracted using the analysis of The A-dependence of 
both pion and nucleon semi-inclusive diffraction is reproduced well by our calculation. A 
priori, the A-dependence for a given semi-exclusive channel could be different from that of 
(Jdiff- However, if fluctuations near the average value of a dominate then, as we discuss 
below, the A-dependence given by eq.2 is sensitive mainly to the value of a to t, see eq.4. Our 
numerical results (using P(cr)) for <7<&//(A) are reasonably close to the calculation of the 



preQCD model of [TS|] because similar values of u a are used. But Ref. |0| showed that their 
calculation agrees with the only experimental data available for the inelastic diffractive total 
cross section, obtained using emulsion targets. Hence we also agree with these data. 

The similarity of the results for different models of P{cr) suggests that fluctuations of 
a near the average give the major contributions. This similarity is even more apparent in 
Fig. 2, where we plot adiff(A)/u a . To test this idea we compute an approximate diffractive 
cross section cr^fff by using a Taylor series about a = (a) in the integrals / daP(a)f(a). 
The result for c^fJJ is 



= I d 2 BT 2 {B)e -(«)T(B)^ (4) 



a diff j 

which is shown as the dashed line in Fig. 2. We see that the approximation is quantitatively 
accurate for A < 50 and qualitatively good for all values of A. This is because for realistic 
P(cr) the dominant contribution to the inelastic diffraction cross section arises from impact 
parameters B near the nuclear surface where (a)T(B) is small. As a result, the second 
cumulant (dispersion of the cross section) dominates the diffractive cross section. This 
shows that the A dependence is mainly determined by the value of (a). The deviation 
at large A must be due to configurations further from those of average cross section; in 
particular the ones of relatively small a (and therefore small size). At the same time our 



numerical analysis shows that the series given by the sum of cumulants of the cross sections 
is badly convergent. This is similar to the poor convergence of the standard Glauber series 
for a tot {hA) which has an A-dependence of A 2 ! 3 while the first term ~ A. 

In the limit that A becomes infinite, configurations of small size can be expected to dom- 
inate since the nucleus acts as a black disk for all other configurations. Indeed Ref. p5]f2T 



suggested that the effects of such small-sized configurations would dominate the pion-nucleus 
inelastic diffractive cross section. This early result is inherent in our eq. (2). The integration 
over small values of a in eq. (2) gives a result oc 1/T(B) for values of B within the nucleus. 
The integration over d 2 B leads to <J di ff oc Ra oc A 1 / 3 . Our realistic functions -P(cr) do not 
vanish at a = 0, so we search for the dominance of small-sized configurations simply by 
increasing A. Numerical evaluations of eq.(2) show that the behavior is close to A 0,33 for 
fantastic values of A greater than about 10000. An additional result indicating that small- 
sized configurations play a small role is obtained by simply cutting off the integrals over a 
at a maximum value of cr max = 5mb. This contribution varies from 2% to 5% as A increases 
from about 12 to 200. This can be considered as an upper limit on the pQCD contribution 



suggested in pi], 25 



The present result does not mean that small-sized configurations are impossible to find. 
Reactions in which the pion diffractively dissociates into two jets of high relative trans- 
verse momentum select those configurations p6] and leads to an A 2 variation of the forward 
diffractive cross section. One could also look for the transition to this A 2 regime by consid- 
ering production of states where the value of M^// is mainly determined by the transverse 
momenta of produced hadrons. 

The same formalism and -P(cx) used to obtain <Jdiff{A) for pion and proton projectiles 
also allow us to compute the A-dependence of the zero angle differential cross section for 
coherent nuclear diffractive dissociation as well as the total cross section. The diffractive 
angular distribution is related to the square of the scattering amplitude M. ~ / d 2 Be l ^' B (h \ 
F(a, B) | X). Squaring M. and summing over the diffractively produced states X yields the 
angular distribution. The result is 
da, i: - 



dt 



l, ff (0°) = 7T / daP{a)f 2 {a) - tt / daP{a)f{cx 



(5) 



where f(a) = J BdB ( 1 — e 2 T ( B n. Numerical evaluation yields the result that dadx}} (0°) 



varies approximately as A 1 ' 24 for the pion (see Fig 3) and A for the proton. The total cross 
section a tat (A) is given by the expression 

atot(A) = 2 J daP(a) J d 2 B <h \ F(a, B)\h> . (6) 



Color fluctuations (also known in this inelastic shadowing have fairly small 

effects on total cross sections. Numerical evaluation shows that the results of using the 
above equation are similar to those of the more detailed calculations of Ref. ||27|| . 



We also need to discuss the energy region for which our calculations are valid. If Eh 
is much less than about 100 GeV, the cross section adiff(A) may increase significantly 
with energy due to the effects of the nuclear longitudinal form factor Fa- Such effects are 
omitted here, but are important in computing the inelastic shadowing correction to the total 
nuclear cross section jETl if Eh < lOOGeV^. In diffractive dissociation of a projectile h into 
a state of mass Mj the minimal longitudinal momentum transferred to the target is 
• F° r sma U Ql the nuclear form factor can be described using the parameterization 
Fa{q 2 ) — exp[— q 2 R\/3}. Thus if one uses P(<r) for energies that vary as A 33 one can 
effectively use the same longitudinal form factor for different nuclear targets. This effect 
is a small correction if Eh is greater than about several hundred GeV. But this is just the 
energy range where diffractive data and data on inelastic shadowing corrections to the total 
cross section of hadron-deuteron scattering are available. A specific energy dependence of 
P(cr) expected at much higher energies is discussed in ref. [|l"3|j . 

Our results show that the effects of color fluctuations are important for coherent nuclear 
diffractive scattering. This has a broader significance because of its relevance in heavy ion 
collisions. In particular, the realistic P(cr), used here to describe diffractive processes, corre- 
spond to a significantly larger probability for multiple scattering processes to occur than the 
usual Glauber approximation. Previous work has shown that color fluctuations leads to 
significant fluctuations of transverse energy in nucleus-nucleus collisions in agreement with 
current data. Moreover the work of Ref. |3j shows that the probability for percolation phase 

f°° dyy n P(y) 

transitions depends strongly on the quantities p n {&) = — t^t\ • These give the probabil- 
ity that n nucleon-nucleon inelastic collisions occur with a cross section larger than a. The 
differences between the different versions of -P(cr) are large and influence the predictions of 
whether or not a percolation phase transition could occur in heavy ion collisions. Measuring 



(Jdiff(A) for proton beams would strongly constrain P(a). 

Color fluctuations have an important intrinsic interest through their close relation to 
QCD. New measurements of coherent nuclear diffraction could determine finer details of 
P(ct) and therefore have a wide impact for studies of heavy ion collisions. We stress that 
new experiments are possible during the next fixed target run at FNAL for a wide range of 
projectiles. 
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Figure Captions 



Figure la. adiff{A), pion beam, for the different P(cr) of \T5\. The data are from Ref. EB 



Figure lb. adiff(A) for a proton beam. The data are from Ref. [24 



Figure 2. adiff(A)/u) a (pion beam). Solid curves from eq.(2). Dashed curve- from eq.(4). 
Figure 3. n forward coherent nuclear diffraction obtained using different P(cr) of JT5[ . 
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